Wave Packet Spreading: Temperature and Squeezing Effects with 
Applications to Quantum Measurement and Decoherence 
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cn ! Abstract 

A localized free particle is represented by a wave packet and its motion is 
l/-) ■ discussed in most quantum mechanics textbooks. Implicit in these discussions 

o . 

is the assumption of zero temperature. We discuss how the effects of finite 
. temperature and squeezing can be incorporated in an elementary manner. 



The results show how the introduction of simple tools and ideas can bring 



the reader into contact with topics at the frontiers of research in quantum 

& . 

mechanics. We discuss the standard quantum limit, which is of interest in the 



measurement of small forces, and decoherence of a mixed ( "Schrodinger cat") 
state, which has implications for current research in quantum computation, 
entanglement, and the quantum-classical interface. 

I. INTRODUCTION 

The uncertainty principle in quantum mechanics implies that the position and momen- 
tum of a particle cannot be determined simultaneously with arbitrary precision. More 
explicitly 

AxAp > ~, (1) 
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where (Ax) 2 — ((x — (x)) 2 ) and (Ap) 2 — ((p — (p)f) are the variance of the position and 
momentum of the particle, respectively. For a noninteracting particle the expected value at 
time t of an operator O is given by 

/oo 
dxil>*(x,t)Oil>(x,t), (2) 
-oo 

where the wave function ip(x,t) is the solution of the time- dependent Schrodinger equation, 

ih^ = H^, (3) 

with H the particle Hamiltonian. For the more general case, such as a particle interacting 
with a heat bath, the expected value would be given by 2 

(0(t)) = Tr{p(t)0}, (4) 

where the density matrix pit) is the solution of the von Neuman equation, 

ih% = [H, P ], (5) 

with H now the Hamiltonian for the entire system. 

For a free particle, the stationary solutions of the Schrodinger equation are plane waves, 
for which the particle may be found with equal probability anywhere in space, that is, 
Ax = oo, Ap = 0. However, one is often interested in describing a localized particle, which 
can be achieved by constructing a wave packet corresponding to a superposition of plane 
waves. Such a packet is necessarily not stationary and will spread (or shrink) in time. Wave 
packet spreading is of fundamental interest, appears in many contexts, and is discussed in 
introductory 1 and advanced 3 ' 4 quantum mechanics textbooks. 

We begin our discussion in Sec. II, where we describe the motion of an arbitrary free- 
particle wave packet. This description is more or less standard, the main result being a 
general expression for the width at time t in term of the initial data. This expression does 
not, however, take into account the uncertainty principle; to do so one must evaluate the 
initial data using Eqs. (2) or (4). We do this in Sec. Ill, where we begin with a brief derivation 



of the uncertainty principle. There we introduce annihilation and creation operators that 
are analogous to those appearing in discussions of the harmonic oscillator, 1,3 ' 4 but which 
now apply to an arbitrary state. An immediate consequence is a simple construction of the 
most general minimum uncertainty wave packet. We conclude Sec. Ill with an expression 
for the spreading of an initial minimum uncertainty state. 

In Sec. IV we consider the effect of finite temperature on wave packet spreading. An 
initial wave packet at finite temperature is in a mixed state: there is no corresponding wave 
function, and the state is described by a density matrix. Nevertheless, we can calculate the 
temperature effect by forming the observed quantities (the probability distribution or the 
expected values) with an initial wave function and then averaging over a Maxwell distribution 
of initial velocities. The result is an additional spreading that dominates when the thermal 
de Broglie wavelength is small compared to the initial width. 

In Sec. V we consider spreading for squeezed states. Squeezing is generally discussed in 
the context of the harmonic oscillator, but here we discuss squeezing of a Gaussian wave 
packet. An interesting result is that, while for sufficiently long times the root-mean square 
(rms) width of such a wave packet increases linearly with time, for short times it can even 
shrink. Finally, in Sec. VI we consider two topical applications: the standard quantum limit 
arising in connection with the measurement of small forces and decoherence. In particular, 
we give a new and simple demonstration of how to circumvent the standard quantum limit. 
This limit is of interest not only for gravitational wave detection, but for any application 
where the question of the accuracy of successive measurements arises. With regard to the 
decoherence problem, the temperature effect on wave packet spreading is an essential feature. 

II. FREE PARTICLE 

We begin by reminding ourselves that from either Eq. (2) or (4), we can show that the 
rate of change of the expected value of an operator O, with no explicit time dependence, is 
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given by: 



it^- = ([0,H\). (6) 



For a free particle with H = p 2 /2m, we can use the canonical commutation relation, 

[x,p]=ih, (7) 

to show that 

d{x) = (p) dM_ = Q , g . 
dt m ' dt 

These are the classical equations of motion (Ehrenfest theorem). Therefore, we have 

<*(*)> = (x(0)> + ^>f, (p(t)> = (p(0)> . (9) 
m 

In the same way, we can show that 

d (x 2 ) (xp + px) d (xp + px) 2 (p 2 ) d (p 2 ) 



dt m dt m ' dt 

and therefore, 



0, (10) 



We can write the above results in terms of the variances as 

AAt) = a^(o) + ±AAo)S + Wp(°) + p(°M°»-2W)W», (12) 

Hence, for sufficiently long times, Ax(t), the rms width of the wave packet at time t, increases 
linearly with time. However, it is possible for a wave packet to shrink for a time, as we shall 
discuss in Sec. V. 

In the formal solution of the equations of mean motion, it appears that the initial data, 
Ax 2 (0), Ap 2 (0) and (x(0)p(0) +p(0)x(0)) - 2 (x(0)) (p(0)), could be given arbitrary values. 
We emphasize that this is not so, and the initial expectations must be obtained from the 
initial state by an expression of the form (2) or (4). In particular, the initial variances must 
satisfy the uncertainty principle (1). 



III. MINIMAL (NON-SQUEEZED) GAUSSIAN WAVE PACKET 



To make the discussion as simple as possible, we begin by restricting our discussion to 
states for which (x) = and (p) = 0. We then introduce the operators 

x op f x .op 

a= Ya +l ^ a= Ya~ % ^ (13) 

where a is a real parameter. These operators are formally identical with the annihilation 
and creation operators usually introduced in connection with the harmonic oscillator, 1 ' 3,4 
but here they apply to an arbitrary state (pure or mixed) without reference to an external 
potential. Next, we form the necessarily positive quantity, 

/ + \ Ax 2 o 2 Ap 2 1 , 

where we have used the canonical commutation relation in Eq. (7). We seek the minimum 
of this quantity with respect to variations of a 2 , which occurs when 



With this value of a 2 , we see that 



AxAp 1 



a'a) = - - > 0, (16) 

which is just the uncertainty principle. The minimum uncertainty state, for which the 
inequality becomes an equality, must be a pure state that corresponds to a wave function <fi 

|2 



satisfying (era] = \\a<f)\\ = 0. That is, must satisfy 

where we have used the familiar realization of the momentum operator: p = The 
solution of the first-order differential equation in Eq. (17) is 

^ = T^W e """^ (18) 

where we have chosen the normalization so that Jf^ dx 4>* (x)(f>(x) = 1. 



To extend this result for nonvanishing (x) and (p), we need only make the replacements 
x — > x — (x) and p — > p — (p) in Eq. (13) and repeat the argument. The result is that the 
most general minimum uncertainty wave packet has the form: 

1 j (x-x ) 2 .mv x\ 

= PF 1 exp I — + *— /' (19) 

where a, xo and vo are real. Thus the minimum uncertainty wave packet is a Gaussian, 
centered at x and moving with velocity vq. 

Suppose we choose the initial state to be a minimum uncertainty state with ip(x, 0) given 
by Eq. (19). Then we find, 

(x(0)) = x , (p(0)) = mv , 

h 2 



x 2 (0)) = x 2 + a 2 , (p 2 (0) =»\ 2 + 



4<r 2 ' 

(x(O)p(O) + p(0)x(0)) = 2mx v . (20) 
With these expressions, Eq. (12) for the mean square width of the wave packet becomes 

A * 2 «=* 2+ yy 2 ' <2i) 

The wave packet will expand so that the mean square width doubles in a time t = 2ma 2 /h. 
During this time, the wave packet will have traveled a distance £ = v^t = 4na 2 /\, where 
A = mv /2nh is the de Broglie wavelength. 

Whereas our derivation might appear to be similar to that found in some textbooks, 
there are important differences in that we have allowed from the beginning the possibility 
of a mixed state and we have shown that the minimum uncertainty state is a Gaussian. 
Furthermore, we have introduced the concept of creation and annihilation operators for 
arbitrary states and, concomitantly (as we shall see in Sec. V), this enables us to consider 
squeezing of arbitrary states (as distinct from just harmonic oscillator states). 
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IV. EFFECT OF TEMPERATURE ON THE SPREADING OF A FREE 

PARTICLE WAVE PACKET 



We proceed by first calculating the wave function for the particle at time t and then 
forming the probability distribution. This procedure is instructive, because it provides 
another method for calculating the result in Eq. (12) for Ax 2 (t). Next, we take into account 
the thermal distribution of initial velocities 

Consider the solution of the free-particle Schrodinger equation with a given initial state. 
The general solution is 4 

/ m r°° , , f mix — x') 2 ) , . , . . . 

= V27rifti Xoo P { 2iht~\ (22) 

We now apply this result to the case of an initial minimum uncertainty wave packet with 
ip(x, 0) of the general form given by Eq. (19). We find, 

1 J (x-x -v t) 2 mv .mv 2 t\ 

The probability distribution is 



P(x;*) = |#r,f)| 



2 



which is a Gaussian centered at the mean position of the particle at time t with variance 
given by Eq. (21). 

These results are standard quantum mechanics. Next we consider an ensemble of particles 
in thermal equilibrium, but so weakly coupled to a heat bath that we can neglect dissipation 
in the equation of motion. Each particle has a wave function of the form (19), with a Maxwell 
distribution of initial velocities. (Note that the wave functions differ only by the phase 
factor exp[imvox/h] and that the distribution in initial velocities implies a corresponding 
distribution of the phase.) We obtain the corresponding probability distribution by averaging 
the distribution (24) over a thermal distribution of initial velocities. The result is 
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I 

: -p(-fe^}, (25, 



^/2ttA4(0 I 2A4(t) 



in which 



Ax 2 T (t) = Ax 2 (t) + — t 2 
m 

= a 2 + f-4r + — 1 **• (26) 

Here we have introduced the subscript T to emphasize that the probability distribution 
corresponds to finite temperature. Thus, there is an additional spreading which is just that 
due to a Maxwell distribution of particle velocities. As we shall see, in Sec. VI B, it is this 
extra term in the spreading of the wave packet which is the origin of decoherence. 

We could just as well obtain Eq. (26) by averaging the expressions in Eq. (20) for the 
initial moments. The result is 

(x(0)) T = x , (p(0)) T = 0, 

2 /J2'~ sX • — 



x 2 (0))^ = x 2 + a 2 , (p 2 (0))^ = mkT + 



t 1 ' \ r v Vt 1 4a 2 ' 

(i(0M0)+K0)x(0)) T = 0. (27) 

If we substitute these results into Eq. (12) for Ax, we find the result (26). 

The center of the distribution Pt(x, t) in Eq. (25) does not move because the mean initial 
velocity is zero. The variance in Eq. (26) is the sum of three terms: the initial variance a 2 , 
the uncertainty principle spreading (ht/2ma) 2 , and the thermal spreading ^t 2 . The ratio 
of the last two is 

™kTa 2 _ a 2 

~ 4A*' (28) 

where A = h/VmkT is the mean de Broglie wavelength. Therefore, the thermal spreading 
will dominate when the mean de Broglie wavelength is small compared to the initial width 
of the packet. 
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Note that the initial thermal state we have described is what is called a mixed state. For 
such a state there is no single wave function with which one can form observables such as 
the probability distribution in Eq. (24); rather the state is described by a density matrix. 2 
We have avoided introducing the density matrix to keep the discussion simple, but for those 
who would like to see it, the matrix is given by 



where is given by Eq. (19). We also note that we have neglected dissipation during the 
time development. Of course, in order to come to thermal equilibrium, a particle must be 
coupled to a heat bath and there must be a corresponding dissipation. The strength of this 
coupling would be measured by a typical decay rate 7. If the coupling is weak, we must 
wait a long time of order 7" 1 for the system to come to equilibrium, but the equilibrium 
state will be independent of dissipation. The situation is like that for an ideal gas: collisions 
are necessary to bring about an approach to equilibrium, but do not appear in the equation 
of state nor in the velocity distribution. On the other hand, the effect of dissipation on 
the time development of the initial state can be neglected only for times short compared to 
7 _1 . Our simple expression (26) for wave packet spreading is therefore valid only for such 
short times (yt <C 1) where the motion is that of a free particle; this short-time behavior 
is exactly what is relevant for the calculation of decoherence times, as we shall discuss in 
detail in Sec. VI B. 

V. EFFECTS DUE TO SQUEEZING 

In general, a squeezed state is defined as one in which the uncertainty of one variable 
is reduced at the expense of an increase in its conjugate variable. 5 If the original state is 
a minimum uncertainty state, then the squeezed state may also be a minimum uncertainty 



<*|p(0)|x'> =0(x)0*(x') 




(29) 
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state (the so-called ideal squeezed state, such as a coherent state), but, more generally it 
is not. We start with the simple case of a minimum uncertainty state which is squeezed 
so that the uncertainty in x remains unchanged but the uncertainty in p increases. The 
corresponding squeezed state is also Gaussian, 

1 (l-iC)x 2 

Mx) = (W)I7I e "^ (30) 

This squeezed state can be represented as the result of a unitary operation on the minimum 
uncertainty state (18), 

4> c {x) = e ^ c(a+at)2 0(:r)- (31) 
If we use the general formula (2) with ip(x, 0) = (f>c(x), we find 

<x(0)> = 0, (p(0)> = 0, 

* 2 (o)W, (p\*)) = h2{l + C2) 



Aa 2 ' 

(x(0)p(0) + p(0)x(0)) = hC. (32) 
The squeezed state is therefore not a minimum uncertainty state, because AxAp = 

IVTTc 2 > f. 

With these expressions, Eq. (12) for the variance of the wave packet becomes 

AAt) = a >( l + ^L)\(J!L)\ ,33, 

W V 2a 2 m) \2maj V ; 

If C < 0, the wave packet first contracts, then expands; for very long times the wave packet 
always expands. 

Finally, it is clear from the analysis given in Sec. IV that the thermal contribution to the 
spreading is the same for both the squeezed and unsqueezed states. 



10 



VI. APPLICATIONS 



A. Standard Quantum Limit 



Accurate measurements of the position of a free mass is a subject of much current interest, 
in particular in the context of gravitational-wave detection where questions have arisen 
as to whether there are fundamental quantum mechanical limits on detection sensitivity. 
The point is that the act of measurement introduces spreading which affects subsequent 
measurements. In this context, Braginsky and Vorontsov 6 have argued that in two successive 
measurements of the position a; of a free mass m made at a fixed time interval t, there is an 
uncertainty Ax(t) in the result of the second measurement satisfying 



The inequality (34) is called the standard quantum limit. On the other hand, Yuen, 7 while 
agreeing that this is the correct result for free masses prepared in coherent states (i. e., 
minimum uncertainty Gaussian states), has argued that the inequality can be violated for 
squeezed states. We now give a simple derivation of the standard quantum limit and Yuen's 
result. 

Consider the squeezed state (30) for which Ax(0) = cr and Ax 2 (t) is given by Eq. (33). 
What is the value of a 2 for which Ax 2 (t) is a minimum for a given tl By calculating the 
derivative with respect to a 2 , we find that the minimum occurs when 




(34) 




(35) 



At the minimum, Ax 2 



(t) has the value 



Ax 2 



(*) 



(VTTc^ + c) 



ht 



(36) 



mm 



rn 



We therefore have in place of (34) the general inequality 




(37) 
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For a minimum uncertainty state, where C = 0, this corresponds to the standard quantum 
limit (34). On the other hand, for C large and negative, the right-hand side of (37) can be 
as small as one likes, as noted by Yuen. 7 

While our explicit result (37) supports Yuen's general conclusion that more sensitive 
detection is possible when the initial measurement results in a squeezed state, a more careful 
examination of the result makes it clear that the may be difficult to achieve in practise. This 
is because the condition (35) puts a restriction on a ( the width of the initial wave packet). 
In fact, a large negative C implies that a be large, so the initial measurement must have 
large uncertainty. Moreover, we see from (32) and (35) that, when the minimum value of 
Ax(t) is achieved, 



Decoherence refers to the destruction of a quantum interference pattern and is rele- 
vant to the many experiments that depend on achieving and maintaining entangled states. 
Examples of such efforts are in the areas of quantum teleportation, 8 quantum informa- 
tion and computation, 9 ' 10 entangled states, 11 Schrodinger cats, 12 and the quantum-classical 
interface. 13 For an overview of many of the interesting experiments involving decoherence, 
we refer to Refs. 11 and 14. 

Much of the discussion of decoherence 15-18 has been in terms of a particle moving in 
one dimension that is placed in an initial superposition state (Schrodinger "cat" state) 
corresponding to two widely separated wave packets, each of the form (19) but having 
Xq = ±| so that the packages are separated by a distance d. Thus, in an obvious notation 
we write the wave-function of the two-Gaussian state as 




(38) 



Thus the energy needed to produce the state increases with increasing C. 



B. Decoherence 




(39) 
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where ipi and ip 2 are each given by the right-side of Eq. (23), but with x replaced by | and 
— |, respectively, and the normalization constant is 

N = 1 (40) 

Hence 

t) = N 2 (\^\ 2 + IV2I 2 + 2Re{^> 2 }). (41) 

Thus, the probability distribution consists of three contributions, two of which correspond 
to the separate packets, whereas the third is an interference term. For a free particle at rest 
at zero temperature (i. e., in (23) we take v — 0), an elementary calculation leads to the 
result 



P(x,t) 



N 2 ( j (x-jf \ j (x + jf \ 

^AxHT)\ V \ 2A^(t)/ +exp \ 2As*(f)J 

d 2 1 htdx \ 



2 J2 



+ 2eXP | 2Ax 2 (t) 8Ax 2 (t)j C ° S Ama* Ax^t))' (42) 
The first two terms are of the single-Gaussian form given by (24) while the interference 
term is characterized by the cosine factor. The key point to be made here is that this 
interference term persists for all time. More generally, when either temperature or dissipative 
effects are present, one measures the disappearance of the interference term i.e. the loss of 
coherence (decoherence) by defining an attenuation coefficient a(t), which is the ratio of the 
factor multiplying the cosine to twice the geometric mean of the first two terms. From an 
examination of Eq. (42), we see that a(t) = 1, corresponding to the absence of decoherence. 

To take into account the effect of finite temperature, we first form the probability distri- 
bution (41) with ipi and ^ 2 given by the form (23) with x = ±| but keeping the terms with 
Vq. It is a simple matter to see that the resulting probability distribution can be obtained 
from (42) with the replacement x — > x — v t. We then average this probability distribution 
over a thermal distribution of initial velocities, as in (25), to obtain the probability distri- 
bution corresponding to a finite temperature T. After a bit of algebra, the result can be 
written in the form 
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Pr(x; t] = = r p r^4w ) + exp v^m) 

x 2 d 2 htdx \ 

+2a(t) eM ~2AxW) " } C ° S WMW J ' (43) 

where the attenuation coefficient a(t) is given by 

— t 2 d 2 
= ^sJaxW) } 

- exp -^^TW ' (44) 

Once again, we see that a(t) = 1 for T = 0. However, for non-zero T and short times (char- 
acteristic of decoherence time scales), whereas the t dependent terms in the denominator 
are negligible, the t dependent terms in the numerator remain and thus we obtain 



a(t) = e- t2/r *, (45) 



where the decoherence time is 



Td = ~od~' ^ 

and v = \j~kT jm is the mean thermal velocity. This is consistent with the results obtained 
in Ref. 17-19 where we have found that the dominant contribution to decoherence at high 
temperatures (kT >> tvy, where 7 is typical dissipative decay rate), is independent of 
dissipation. However, for very low temperatures T, dissipation plays an important role, in 
which case one must use sophisticated techniques from non-equilibrium quantum statistical 
mechanics. 16,18 

In order to see why decoherence is a short-time phenomenon, consider as an example 
an electron at room temperature (300 K), v — 6.8 x 10 6 cm/s. Hence, if we take d — 1 
cm and a = 0.4A, then using Eq. (46), we obtain = 6.9 x 10~ 24 s, which is orders of 
magnitude smaller than typical 7 _1 values. For this reason it is permissible to take jt « 1 
for calculations involving the calculation of decoherence times and this is why the simple 
derivation outlined above (which is solely within the framework of elementary quantum 
mechanics and equilibrium statistical mechanics) works. 
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VII. CONCLUSIONS 



Wave packet spreading is of fundamental interest in quantum mechanics. By extending 
some of the usual methods, we have been led to many interesting phenomena: tempera- 
ture and squeezing effects on wave packet spreading with applications to topical phenomena 
such as how the standard quantum limit may be circumvented by squeezing and how the 
temperature dependence of the rate of decoherence may be calculated in a simple man- 
ner. Concomitantly we have introduced simple tools and ideas that are at the frontiers of 
cutting-edge research in quantum mechanics, such as Schrodinger cat states, decoherence, 
entanglement and the classical-quantum interface. 
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